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ABSTRACT 

We present the first interferometer map of the H 2 maser emission in 
the active nucleus of NGC4945, which exhibits both starburst and Seyfert 
qualities. Although the declination of the galaxy is —49°, we were able to 
make the observations with the southernmost antennas of the Very Long 
Baseline Array. Strong maser emission is present in three velocity ranges, 
one near the systemic velocity and two shifted roughly symmetrically by 
±(100 — 150) km s _1 . This is the first detection of highly blue-shifted water 
emission in NGC4945. We also report a marginal detection of previously 
unreported red-shifted emission at a velocity of ~ +200 km s" 1 with respect to 
systemic. From fringe rate analysis we determined the position of the maser to 
be a 1950 = 13 A 02 m 32f28 ± 0f02; S 1950 = -49°12'01"9 ± O'/l. The uncertainties 
in earlier estimates are at least several arcseconds. The maser lies within 2" 
(36 pc at a distance of 3.7 Mpc) of the peaks in 1.4 GHz continuum and 1.6 
/iin emission from the nucleus. The mappable maser emission is distributed 
roughly linearly over ps 40 milliarcseconds (0.7 pc) at a position angle of ~ 45°, 
which is close to the 43 ± 2° position angle of the galactic disk. The red and 
blue-shifted emission symmetrically stradle the systemic emission on the sky, 
which suggests material in edge-on circular motion around a central object. The 
position- velocity structure indicates a binding mass of ~ 1 x 10 6 M Q , within a 
volume of radius ~ 0.3 pc, which in combination with estimates of the AGN 
luminosity, implies that the central engine radiates on the order of 10% of its 
Eddington luminosity. 

Subject headings: galaxies: individual, NGC4945 -- galaxies: kinematics and 
dynamics — galaxies: nuclei — masers 
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1. Introduction 

The edge-on spiral galaxy NGC4945 contains the first H 2 maser discovered in a 
galactic nucleus (|dos Sautos fc Lepine 1979] ; [Batchelor, Jauncey, fc Whiteoak 1982|) . Within 



the nucleus there is a heavily absorbed hard X-ray component (> 10 keV), with a measured 
column density of nu = 5 x 10 24 cm -2 ( [Iwasawa et al., 1993|) , and true X-ray luminosity 



(2-20 keV) of 9 x 10 41 erg s" 1 . The hard component varies on timescales of hours, implying 
a source size of ~ 10 15 cm. The nucleus is very prominent at infrared wavelengths, with a 
far infrared luminosity Lfir = 2 x 10 43 erg s _1 from the inner 200 pc ( [Brock et al., 198lf ), 
at an assumed distance of 3.7 Mpc flWhiteoak fc Wilson 1990|) . The nucleus also shows 
a~lx6 kpc ionization cone along the galactic minor axis (|Nakai 1989Q and emission 
knots with LINER-like spectra ( |Moorwood et al., 1995| ). Optical line splittings of up to 



600 km s (|Hcckman, Armus, fc Miley 1990| , [Moorwood et al., 1995|) indicate that the cone 



is evacuated by a nuclear wind. While X-ray characteristics suggest that the active galactic 
nucleus (AGN) is driven by an obscured Seyfert nucleus, the high IR luminosity and wind 
indicate the presence of a nuclear starburst. 

Spectra of the maser consistently show emission (~ 1 Jy) within about ±20 km s~ of 
the systemic velocity and stronger emission (~ 5 — 10 Jy) red-shifted by ~ 100 — 150 km s" 1 
(e.g., Whiteoak fc Gardner 198~6] , Figure 1). Previously, no corresponding emission on the 



blue side of the systemic velocity has been detected (e.g., [Nakai et al., 1995] ). We adopt 
a heliocentric systemic velocity of 561 ± 4 km s _1 , based on CO observations ( pahlem~et| 
al., 1993| ). (Velocities are defined with respect to the standard radio astronomical definition 



of Doppler shift throughout.) The importance of nuclear H2O masers to the study of 
AGN has been widely appreciated after recent observations of the maser in the nucleus of 
NGC 4258 ([Miyoshi et al., 1995|) , where the masers were found to trace a thin molecular 
Keplerian disk of 0.13 pc inner radius around a central mass of 3.6 x 10 M . In this Letter, 
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we report the first sub-arcsecond estimate of the maser position in NGC4945, the detection 
of blue-shifted maser features, and the first interferometric map of the distribution of maser 
emission. 

2. Observations & Results 

We observed the maser with the Very Long Baseline Array of the NRAO 1 for 4 hours 

1 The National Radio Astronomy Observatory is operated by Associated Universities, Inc, 
under cooperative agreement with the National Science Foundation. 

on 1995 January 22. We used the southernmost stations, at Mauna Kea, HI (MK), Saint 
Croix, VI (SC), and Kitt Peak, AZ (KP). At KP the source reached a maximum elevation 
of 8.6°, while it was 11° and rising at MK. We obtained useful data for only 2.5 hours on 
the MK-KP baseline. Observing conditions at Saint Croix were relatively poor and no 
useful data were obtained from it. Tipping scans showed that the opacity at KP, MK, and 
SC were 0.07, 0.05, and 0.3 at zenith, respectively. The baseline length was 4470 km and 
the angular resolutions were ~ 0.6 and ~ 1.8 milliarcseconds (mas) in right ascension and 
declination, respectively. The system temperatures for 8° elevation at KP and MK were 
170 and 130 K, respectively. We recorded four contiguous bandpasses covering the range 
405 km s _1 to 821 km s _1 , so as to include the systemic emission, the red-shifted emission, 
and the likely velocity range of counterpart blue-shifted emission. The spectrum is shown 
in Figure 1. 

We calibrated the data with the NRAO AIPS package. The interferometer delays and 
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fringe rates were estimated from observations of the continuum sources 4C39.25, 3C273, 
and 1147 — 382. The latter source was observed at elevations < 20° and was included 
so that we could estimate the magnitude of elevation dependent errors in the correlator 
model and post-correlation analysis. The delay calibration is accurate to about 3 ns overall, 
after implementing the U. S. Naval Observatory station positions and earth orientation 
parameters (Eubanks, private communication). 

We analyzed the fringe rate of the maser emission as a function of time to estimate 
the position of the maser, taking into account the effect of the atmosphere. We fit a 
sinusoid to the measured fringe rates for emission between 664 and 719 km s _1 (with 
uniform weighting) to obtain position offsets of +2.7" and +8.2" from the a priori position 
^1950 = 13 ft 02 m 32f000, 5i95o = — 49°12'10"05, with respect to a reference frame consisting of 
the sources 4C39.25, 1147 — 382, and 3C273B. We conservatively estimate the uncertainties 
in position based on the residual fringe rates of the maser and calibrator sources. The new 
maser position is 

a 1950 = 13 /l 02 m 32f28 ± Of 02, 5 1950 = -49°12'01'"9 ± 0"l; 
a 2000 = 13 /l 05 m 27f48 ± Of 02, S 200 o = -49°28'05"4 ± 0"1. 

To map the brightness distribution of the maser, we phase referenced the data to the 
strong spectral feature at 699 ± 0.8 km s~~ . The phases of this feature were subtracted from 
the phases of other spectral features, which corrected for the effects of the atmosphere and 
calibration uncertainties. We assumed a point source model and fit a sinusoid (of 24 hour 
period) to the difference phase, estimating relative feature positions and formal errors (e.g., 
Thompson, Moran, fc Swenson 1986| ). Sample phase plots are shown in Figure 2. The data 
were averaged in frequency across the line profile for each feature and averaged in time for 
between 300 s and 640 s so that phase ambiguities could be resolved prior to fitting. For 
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a 1 Jy feature that is 1 km s _1 wide (half power width) the signal-to-noise ratio in 300 s is 
about 4. 

The structure of the maser is remarkably linear, with emission distributed over « 40 
mas between 414 and 719 km s _1 (Figure 3). Deviations from linearity are antisymmetric 
with respect to the systemic emission, which gives the apprearance of a shallow "S" on the 
sky. The centroids of the blue and red emission are separated by ~ 30 mas or 0.55 pc, at a 
position angle of s=s 45°, and roughly symmetrically bracket the systemic emission. The wide 
separations among the systemic, red, and blue emission correspond to many interferometer 
fringes and are well determined. Maser features that are detected only marginally (i.e., with 
a rms of residuals about the fitted model of > 1 rad) are not shown in the map, although 
they do not appreciably alter the source structure. These marginal features include one 
at ~ 779 km s -1 , which is more red-shifted than any previously reported maser emission 
in NGC4945. In Nonetheless, not all of the detected maser features could be mapped, 
particularly near the systemic and blue-shifted velocities, because a significant fraction are 
not well approximated by a point source model (i.e., the phases could not be fit by a diurnal 
sinusoid) . 

The interferometer calibration was sufficiently good that uncertainty in the astrometric 
position of the maser is the largest source of quantifiable systematic error in estimating 
maser feature positions. The error is <^ 90 yuas, which is computed for an astrometric 
uncertainty of 100 mas and for a feature 200 km s _1 from the reference, and is proportional 
to both quantities. (Residual atmospheric delay contributes a moderate systematic error 
because of the low elevations at which we observed the maser; it is smaller than the error 
from the astro metry.) A second source of systematic error arises from the possibility that 
emission, for each feature, may not be well described by a point source, which is difficult 
to quantify. However, the data for each mapped maser feature is well fit by a point source 
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model over the limited range of hour angle available. The systematic error in position for 
each maser feature is less than the statistical error, which is dominated by the limited range 
of hour angle over which the phases were measured. This limitation leads to 1) a large 
correlation between right ascension and declination parameters and 2) an uncertainty as 
to in which lobe of the interferometer fringe pattern the maser lies. Uncertainties in the 
estimated positions take into account possible lobe ambiguities. 

3. Discussion 

The new astrometric position of the maser clearly shows that it is associated with the 
center of the galactic nucleus (Figure 5). The previous position of Batchelor et al., (1982), 
obtained with the Parkes antenna of the CSIRO, was substantially more uncertain because 
of the beamwidth of the antenna. The maser lies within la of the infrared peak at 1.6 ^m 
and overlaps the radio continuum source, which Whiteoak & Wilson (1990) estimate is 
6" x 2" in size at 6 GHz (see Figure 5). 

We cautiously interpret the data for the maser in NGC 4945 in terms of the paradigm 
established for NGC 4258. The extensive data on NGC 4258, with its clear Keplerian 
signature, show that the masers define a nearly edge on thin annular disk. The data for 
NGC4945 are more sketchy but they share the characteristic of systemic emission bracketted 
roughly symmetrically in velocity and position by high-velocity emission. Moreover, the 
angular distributions of blue-shifted emission is roughly linear and the velocities decline with 
angular distance from the systemic emission. However, the position errors are significant, 
the red-shifted emission does not appear to mirror this structure, the decline in velocity for 
the blue-shifted emission is faster than what is characteristic for a Keplerian rotation curve, 
and the non-point source emission that seems to exist at some systemic and blue-shifted 
velocities could not be mapped (Figure 4). Nonetheless, there is a suggestion that the 



structures of NGC4945 and NGC4258 may be similar, and we proceed on that assumption. 

In the context of an edge-on disk, the observations suggest a binding mass (for circular 
orbits) of M — l.iy 2 £>Mpc^mas esc 2 i esc 3 9, where M is the enclosed mass in units of M Q , V 
is the orbital velocity in units of km s _1 , D Mpc is the distance in Mpc, 6 mas is the radius in 
mas, i is inclination, and 9 is the azimuthal angle of the red-shifted masers in the disk plane 
(90° is the plane of the sky). The binding mass isM» 1.4 x 10 6 M Q for V = 150 km s _1 , 
b = 15 mas or 0.3 pc, and i = 9 = 90° (Figure 4). The implied mass density is ~ 1.6 x 10 7 
M pc~ 3 . The position angle of the disk, ss 45°, is comparable to the position angle 
of the larger galactic disk, 43 ± 2°, and the senses of rotation of the two disks are the 
same. If the observed 9 x 10 41 erg s _1 X-ray luminosity ( Iwasawa et al., 1993 ) implies a 



bolometric luminosity on the order of 10 43 erg s _1 , then the central engine radiates 10% 
of its Eddington luminosity. We note that the shallow "S" the mappable maser emission 
traces may be consistent with a warp in the disk, which could be caused by radiation-driven 
torques that result from oblique illumination by the central engine QMaloney, Begelman, fc 
Pringle 1996|) , though this is not a definitive explanation. 



The interpretation of the observed maser structure is not unique and radically different 
interpretations might be viable. The emission could arise from a complex of star forming 
regions superposed on or within the nucleus, or from jet-induced shocks in the AGN. The 
linear extent of the maser is comparable to galactic star formation complexes (e.g., W49N), 
although the velocity difference between the systemic and red-shifted features is large for 
adjacent star forming regions. Maser emission can arise in warm gas behind dissociative 
shocks flKlitzur, Hollcnbach, fc McKee 1989| ) and nondissociative shocks ( |Kaufman &] 



Neufeld 1996| ), although the maser spectrum mandates a complex line-of-sight velocity 



structure for the shocks. 

Interferometric synthesis observations will undoubtedly clarify the dynamical 
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characteristics of the maser region, providing a comprehensive map of the maser emission 
within 200 km s _1 of the systemic velocity. The dual nature of the nucleus, as both 
a starburst and a very heavily obscured Seyfert object, makes this water maser source 
particularly important since it may be used to probe the dynamics and structure of the 
nucleus on sub-parsec scales. 

We thank P. Diamond for his assistance with aspects of the post-correlation processing, 
and S. Ellingsen for assistance in obtaining a spectrum of the maser at Parkes. 
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Fig. 1- Spectrum of the water maser emission from the nucleus of NGC4945 made on 1995 
October 13 at the Parkes radiotelescope of the CSIRO. The maser is variable on time scales 
of months but the spectrum shown is typical. The integration time is about 120 minutes 
(on and off source). The adopted heliocentric systemic velocity of the galaxy, 561 km s _1 , 
is indicated by the vertical bar. The dashed lines indicate the extremes of velocity that 
we observed with the VLBA. The solid bars above the velocity axis indicate the range of 
velocity within which maser emission is detected. The absolute gain calibration is accurate 
to ~ 30 %. Velocities with respect to the Local Standard of Rest may be obtained from the 
relation Vlsr = Vjiei — 4.6 km s _1 . 

Fig. 2- Sample plots of phase vs time for (top) the systemic maser feature at 557 km s _1 and 
(bottom) the red-shifted feature at 665 km s _1 . The phases sample a 24 hour sinusoid and 
are relative to the phase of the reference feature at 699 km s _1 . The right ascension and 
declination offsets of the maser features, with respect to the reference, were determined 
from these curves. The hour angle of the source with respect to the baseline is —8.5 h at 12 
UT. 

Fig. 3- Map of the maser brightness distribution. The gray-scale shows velocity (414- 
719 km s _1 ) with lighter shading indicative of more blue-shifted emission. The error bars 
reflect measurement noise at the 3<r level. The position angle of the galactic disk is 43 ± 2° 
( pahlem et al., 1993| ). Marginal detections, for which the rms of the phase residuals from 
the model are > 1 rad, are also not shown. The arrow indicated the maser feature closest 
to the systemic velocity of the galaxy. 

Fig. 4- Position- velocity diagram for the maser emission. The impact parameter is measured 
with respect to the maser feature closest to the systemic velocity of the galaxy. The dashed 
line indicated the systemic velocity. The error bars represent position uncertainties at the 
3a level. 
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Fig. 5- Comparison of various position determinations for the nucleus of NGC4945. 
These include the new maser position reported here and the a priori position obtained 
by Batchelor et al., (1982) with the Parkes antenna, the interferometric position of H I 
absorption flAbles et al., 1987|) , the synthesis positions at 408 MHz ([Large et al., 1981]) , 1.4 



GHz ( |Ables et al, 1987|) , and 6.0 GHz flWhiteoak fc Wilson 1990]) , and the position of the 



1.6 jum peak (Whiteoak & Gardner, unpublished). The sizes of the symbols correspond to 
the la uncertainties in position. The deconvolved angular size of the continuum source 
at 6 GHz is about 6" x 2" ( |Whiteoak fc Wilson 1990|) at a position angle of ~ 45° and is 



shown by the dashed ellipse. The position angle of the galactic disk is 43 ± 2° ( PahlenTet 
al., 1993|) . 
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Fig. 1. — Greenhill et al. 
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Fig. 2. — Greenhill et al. 
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Fig. 3. — Greenhill et al. 
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